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ABSTRACT 



Context. A derivation of a generalized Ve-law for nonthermal collisional rates of excitation by charged perturbers is presented. 
Aims. Aim of this paper is to find a more general analytical expression for a surface value of the source function which can be used 
as an addtional tool for verification of the non-LTE radiative transfer codes. 

Methods. Under the impact approximation hypothesis, static, one-dimensional, plane-parallel atmosphere, constant magnetic field of 
arbitrary strength and direction, two-level atom model with unpolarized lower level and stimulated emission neglected, we introduce 
the unphysical terms into the equations of statistical equilibrium and solve the appropriate non-LTE integral equations. 
Results. We derive a new analytical condition for the surface values of the source function components expressed in the basis of 
irreducible spherical tensors. 
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1. Introduction 



In the series of papers offLandi Pegi'Innocenti et al.l dl991alfbh . 
lLandi Pegi'Innocenti & Bommiei ~ ?1994l) (from now on refer- 
enced as iPaper J) , the general formalism of resonance line po- 
larization scattering for a two-level atom has been developed. 
The non-LTE problem of the 2nd kind for an arbitrary magnetic 
field, three-dimensional geometry of the medium and arbitrary 
irradiation by external sources has been discussed. The effect of 
inelastic collisions with charged perturbers has been considered 
for the particular case of a relative Maxwellian velocity distribu- 
tion. 

IPaper II analysed the analytical properties of the solutions in 
the particular case of a one-dimensional, semi-infinite, static at- 
mosphere with a constant magnetic field of arbitrary strength 
and direction and assuming zero external irradiation of the at- 
mo spher e. They derived a generalization of the well known -y/e- 
law (e . g. I Avrett & Hummer|[l965l; lMihalas|[l970l: iHubenvll 1 987t) 
for the case of polariz ed radiation and extended the previous re- 
sults of llvanovl(ll990h who studied scattering in a non-magnetic 
regime. 

In most cases of practical interest the polarization degree is 
rather small. The purpose of this paper is to find a new analytical 
solution of the non-LTE problem in unphysical conditions in or- 
der to better verify the accuracy of the polarized radiation trans- 
fer codes. This is done by introduction of an unphysical source 
term in the polarization into the equations of statistical equilib- 
rium. Such a generalization can be useful in testing the accuracy 
of the radiative transfer codes whose purpose is to deal with the 



non-thermal collisional processes (for instance in the impact po- 
larization studies of solar flares). 

Following the approach of the papers quoted above, we adopt 
the formalism of density matrix in the rep resentation of irre- 
ducible tensorial operators (e.g. iFanol Il957l) . We consider the 
lower level with total angular momentum j to be unpolarized. 
This level is completely described by the overall population 
which is set to 1 for normalization reasons. The upper level with 
angular momentum / is described by the multipole components 
of the statistical tensor p^. Coherences between different levels 
j and / are neglected but coherences between Zeeman sublevels 
of level / are in general taken into account. The calculation is 
performed in the Wien limit of line frequency whose assump- 
tion makes it possible to neglect stimulated emission effect, and 
to preserve the linearity of the non-LTE problem. 



2. Equations of statistical equilibrium 

The suitable coordinate system So for atomic state description 
is the one with the z-axis directed along the magnetic field (see 
Figure [TJ. 

Radiative rate contri butions to the evolution of s tatistical op- 
erator Pq are given by (lLandi Pegi'Innocenti! 1985h 



dp K o 



dr 



Send offprint requests to: J. Stepan 



= -iA n TQp K Q -A fj p K Q 



V27TT 



>jj' j -q- 



(1) 



2 



J. Stepan and V. Bommier: Generalized -ye-law 



In this equation Afj (Bjy) is the Einstein coefficient of spon- 
taneous emission (absorption) from level / (J) to level j (/). 
F = 2ngfVL/A/j with g? being the Lande factor of the level 
/ and vl is the Larmor frequency. The transition-dependent nu- 
merical factor wfj has been defined by lLandi DeglTnnocentil 
( 1984) as have the irreducible components of the mean radiation 

tensor J Q . Besides the radiative rates, collisional rates have to 
be considered in the statistical equilibrium, because the source 
of radiation in a semi-infinite atmosphere is the collisional exci- 
tation followed by radiative de-excitation. Thus, the source term 
of the radiative transfer equation originates in the inelastic col- 
lision effect. As the purpose of the present paper is to consider 
unphysical source terms in the non-zero ranks (K, Q) of the ir- 
reducible tensorial operator basis Tq, we will introduce an un- 
physical (K, ^-dependence to the inelastic collisional rates of 
the statistical equilibrium equation below. The purpose here is 
not to thus describe anisotropic collisions, which would require 
a proper formalism that is out of the scope of the present p aper 
(see, for instance, Land i DeglTnnocenti & Landolfil d2004l) for 
a two-level atom, and iDerouichl d2006l) . for polarization trans- 
fer rates in a multi-level atom due to isotropic collisions). We 
introduce as usual the depolarizing rate due to isotropic elastic 
collisions. Thus, the contribution of collisional rates reads 



dp K o 



dt 



(2) 



COLL 



The terms (Cy)g and (Cj)g on the right-hand side of equation 
© are the multipole components of collisional rates of excita- 
tion and relaxation respectively. D (K) is the depolarization rate 
due to elastic collisionsO] 

The radiative and collisional ra tes can be added under the 
impac t approximation hypothesis ( Bom mier & Sahal-BrechoJ 

[Hail) d P K Q /dt = [d P £/df] RAD + [d P K Q idt] coll- Using the equa- 
tions (fl}, (0, and the condition for static atmosphere, dp^/df = 
0, we obtain the equations of statistical equilibrium 

[iA fj rQ+A n + (Cjj) K Q +D^]p K Q 

W (K \-1)Q K 
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By applying the relation between Einstein coefficients for spon- 
taneous emission and absorption, 



2/ + 1 c z 

D JJ ~ o ; , 1 r,,3 J'J 



2j + 1 2hv\ 

and dividing the formula <£3j by Ayj, we obtain the equation 

( l +e K +6 (K) +iTQ)p K 

V27TT c 2 (C jr ) K Q 
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One can introduce the dimensionless parameter of the depolar- 
ization rate 

D (K) 



6 (K) = 



Afj ' 



(6) 



1 This process cannot change a total population of the level. 
Therefore it is always O (0) = 0. We take formally into account only 
the depolarization rate £) <K) to use a formalism coherent with the previ- 
ous papers. A general treatment of physically more relevant transfer of 
multipole components of the upper level is out of scope of this paper. 




Fig. 1. The reference frame Si has its Z-axis oriented vertically 
with respect to the atmosphere, while the z-axis of the reference 
frame Zo is parallel to the direction of magnetic field B. The axes 
X and x lie in the same plane defined by Z-axis and B; the axes 
Y and y are defined to complement the right-handed orthogonal 
coordinate systems. 



and the irreducible tensor which plays the role of generalized 
photon destruction probability 



4 = 



(7) 



If the relation (C^Pq ^ is satisfied we may define the quantity 



B (.KQ) = 2kV l 



2/ + 1 (c 



J'JQ 



(8) 



It is easy to show (see below) that in the particular case of a 
Maxwellian velocity distribution of colliders the relation Z? (00 ^ = 
Bp is satisfied, where B P is the Planck function in the Wien limit 
at given temperature. Using the definition of irred ucible compo- 
nents of the two-level source function (cf. Pa per J) 



s K^hvl 2J+1 ^ 



V27TT 



PQ, 



(9) 



we obtain the statistical equilibrium equations in the compact 
form 



( l +e K +s (K) +iTQ)s K 

= w^(-i)Qj k _ q + 4b^\ 



(10) 
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3. Solution of the Wiener-Hopf equations 

From now on we reduce our analysis to the case of semi-infinite, 
plane-parallel geometry with constant magnetic field along the 
atmosphere. The velocity distribution and volume density of col- 
liders is also constant along the atmosphere but it is in general 
non-thermal. The only position coordinate is the common line 
optical depth t. 

Following the procedure of Paper I a formal solution of ra- 
diative transfer equation is substituted into the definition of ten- 
sor J q, after that we obtain a set of integral Wiener-Hopf equa- 
tions of the 2nd kind, 

{\ + e^ + 5 {K) +irQ)S K Q {T) 

Fq Jo 



We put these results into (TT~5T > to get 



(19) 



Calculation of the limit t — > oo of both sides of the equation 
(fT2l leads to the set of linear algebraic equations for the compo- 
nents of source function tensor in the infinite depth: 



KijWt 

CO 



S j(oa) = bi. 



We can solve these equations and write 
S(oo) = L l b, 



(20) 



(21) 



(11) 



which describe coupling of the tensors Pq(t) at differ- 
ent optical depths via radiation. Several important prop- 
erties of kernels Kko.k'Q'{t,t' ) have been discussed by 
lLandi DeglTnnocenti et al.fi fl990) and in iPaper j Using their in- 
dexing notation one can rewrite the equation ( fTTT i in the short- 
handed form 



where 5 is the formal vector of 5, components, b is the formal 
vector of bi components, and the elements of matrix L are de- 
fined by relation 



{L}ii = a 



oo 



(t)dt. 



(22) 



Establishing a new matrix t — L 1 and substituting (f2TT > into 
( |T9b leads to the generalized form of the -\/e-law 



ai Si(T) = Y f K u (\t - T'\)Sj(T'W + b h (12) 2 aiSi(0)2 = Z W<7- 



(23) 



with 

en = 1 + + 8 (K) + irg, 



bi 



Ji p(KQ) 



(13) 
(14) 



The index i in these expressions runs between the limits 1 and 
N, where N is the number of g-multipoles. In the following 
we br iefly repeat the derivation performed by Frisch & Frischl 
( 1 1975b emphasizing the differences due to presence of bi terms. 

Calculation of the derivative of ( TT2l with respect to r, split- 
ting the integral on the right-hand side into two parts, multiplica- 
tion of the equation by S ,(r), summation over index i, and finally 
integration with respect to t leads to the set of equations 

Si(f)— — At 
At 





Xco 
Kij^SiWAi 



4. Particular solutions 

Setting the special conditions for magnetic field and collisional 
rates, one recovers the less general but more common and ex- 
plicit formulations of the -\/e-law than the one given by ( 1231 . 
In the following sections we will verify this result in the limit- 
ing conditions assumed in recent papers and we will analyse the 
simple examples of non-thermal collisional excitation. 

4.1. Maxwellian velocity distribution of colliders 

In the case of Maxwellian velocity distribution of colliders, re- 
laxation rates of all multipole components are the same: 



(24) 



where , is the usual relaxation rate for collisional deexcitation 



j j 



Xco /-*oo 
A T Si(T) J o At'K u (\t' 



t|) 



I.J 



dSj(T') 
At' 



(15) (C jf ) 



from / to j. For excitation rates one has 
Ci 



Jf'Q 



V27TT 



(25) 



The left-hand side of ( [TBI ) is easily evaluated as 
^ fl ,[5,(^) 2 -5,(0) 2 ] 



(16) 



The first term on the right-hand side of ( fTBI l is evaluated using 
the kernels symmetry Kij(t) = Kji(t) and the equation (fT2b . so 
that we obtain 



where the factor (2/ + l) -1 ^ 2 has been introduced to make a 
connection with the usual collisional rate Cjf of standard unpo- 
larized theory. In this isotropic case, there is no collisional exci- 
tation of higher ranks of density matrix. From the assumption of 
thermodynamic equilibrium one has 



YjSi(0) [aiSi(0)-bi]. 



c*. 

J J 



+ 1 c -hv /k B T 

2 j + 1 



(26) 



(17) 



while the second term equals 

l - at [Si(oo) 2 - Si(0) 2 ] -^i [Si{°°) - S,(P)] . (18) 



where stands for the Boltzmann constant and T for a tem- 
perature of the atmosphere. From ( l24t and ((7]) it is evident that 
€q = e for all possible K and Q, where e is the common photon 
destruction probability. Further 



B {KQ) = B p 6 ko 6qo. 



(27) 
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Substituting the rates ( l24b and d25l l into the formula ( |22] | and 
employing the general identity J^Kn(t)dt = Sa (see lPaper"H) 
together with /?, = we recover form ( l23b the formula (16) of 
the previously cited paper: 



Yjl + e + S (K) + ir0)[S£(O)] 2 = eB\. 



(28) 



KQ 



Assuming that there is zero magnetic field, i.e. F — 0, the 
source function tensor reduces due to symmetry reasons to the 
two non-vanishing components S% and S 1 ? in the reference frame 
Ei . This reference frame is suitable for descriptions of the atomic 
system under these conditions, so that we may identify Eq = Ei, 
with X and Y axes oriented arbitrary in the plane parallel to atmo- 
spheric surface. Further, assuming that there is no depolarization 
of the upper level {6 (K) = 0), we realize from i 



[S»(0)] 2 + [S2(0)] 2 



1 +£ 



-B P 



(29) 



which is the same result derived in different notation bv llvanovl 
( 1990). For simplicity the common alternative to the photon de- 
struction probability has been introduced: e' = e/(l + e). 

If depolarization of the upper level is high enough to destroy 
atomic level polarization (5 (Ar) — > oo for K > 0), or the upper 
level is unpolarizable, the common V^-law for scalar radiation 
is recovered, 



1 + e 



-Bp = Ve'Bp. 



4.2. Anisotropic alignment (de)excitation 

The relation e^ = e is not in general satisfied for all the mul- 
tipoles because the relaxation of the state depends on the 
velocity distribution of colliders. In the following text we will 
neglect the effects of magnetic field. 

Let us assume an example of a relative velocity distribution 
of particles that is axially symmetric with the axis of symmetry 
parallel to the vertical of the atmosphere (so that it is as in the 
former case Eo s Ei) and that the collisional interaction can be 
fully described by only the first two even multipole components 
of this distribution. Thanks to these assumptions the only non- 
vanishing excitation collisional rates are (C//<)2 and (C//)g, the 
relaxation rates (Cy ■){] and (C^) 2 , and for the same reasons the 

only non-zero source function components are and Si. 

An explicit evaluation of the integrals of kernels 
f KicQ,K'Qi(T,T')dT' under these conditions shows that 
the only non-zero ones are giv en by (A5) and (A 12) of 
lLandi DeglTnnocenti et al . ( 1991b). In our notation they read 



f 



^oo,oo(T,r')dr' 



^20,2o(T,r')dr' = — W 2 , 

„( 2 K2 



(31) 



(32) 



with W2 = (wyp 2 . Substituting these results into d23l we see that 

(I + e° )(S° ) 2 + (l+e 2 )(S 2 ) 2 

(e 2 B (20) )2 



e Q Q (B m ) 2 + 



(33) 



To check the validity of polarized radiative transfer codes, 
it is advantageous if one can verify that the transfer of higher 



ranks of the radiation tensor is accurate enough. In the realistic 
scattering polarization models the polarization degree does not 
exceed a few percent so that |Sq(0)| » [5g(0)|. By setting ar- 
bitrary (even unphysical) collisional rates it is possible to verify 
transfer codes in conditions with \S^\ <K \Sq\. 

To privilege transfer in higher ranks of the radiation tensor 
one can artificially suppress the excitation rate (Cjf)^. In the 
extremal case one can set (C/f)o ~ > 0. The easiest way to do 
this is the formal interchange of the role of excitation rates of 
population and alignment, i.e. (C//)2 <-* (C//)q of the original 
Maxwellian velocity distribution: 



(c jf r = 0, (CjjX = 



2 _ ^JJ 



V27TT 



(34) 



(no collisional excitation to upper level population) and the re- 
laxation rates set to the Maxwellian ones 



In this case we have 

g (00) = Qj B (20) = 

and again 



4 = 4 = 



(35) 



(36) 



(37) 



(30) Substituting this into d33l we find out the yfe-law in the form 



2 - [5 2 (0)]2 



^1 - iW 2 (l - 6') 



Bp. 



(38) 



The particular collisional rates d34l > are in fact arbitrary and 
have been chosen to obtain a formula similar to the one of the 
Maxwellian distribution case. 

This relation is useful to test polarized radiative transfer 
codes, because in this unphysical case 5g(0) is the largest term, 
unlike the physical case where the largest term is S jj(0) and S q(0) 
is only a few percent of it. By applying Eq. d3~8l the test is much 
more sensitive to the polarization, and the polarization is better 
tested. We have thus successfully tested a multilevel non-LTE 
radiative transfer code that we are developing, but this code and 
its results are the subjects of a forthcoming paper. 

5. Conclusions 

We have derived a more general formulation of the so-called 
Ve-law of radiation transfer. This analytical condition couples 
the value of source function tensor of a two-level atom with 
other physical properties of the atmosphere. The simplest re- 
sult obtained in conditions of a non-magnetic, isothermal, plane- 
parallel, semi-infinite atmosphere with thermal veloci ty distri- 
butio n of particles and unpol arized a t omic levels (e.g. iMihalas 
119701) has been generalized by Ivanovl ( 119 90) to account for scat- 
tering of polarized radiation and polariz ed upper atomic level. 
Further generalizations done in Paper I, which account for a 
magnetic field of arbitrary strength and direction, has been ex- 
tended in the present paper to account for non-thermal colli- 
sional interactions. It was done by introducing the tensor of the 
photon destruction probability £q and by defining the function 

B (KQ) . 

The resulting formula ( l23l reduces to the cases mentioned 
above if the physical conditions become more symmetric. On 
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the other hand, situations with a high degree of perturbers ve- 
locity distribution anisotropy and especially ones with unphysi- 
cal collisional rates result in a wide range of models which can 
be calculated both numerically and analytically. Thus they offer 
new possibilities for verification of the non-LTE radiation trans- 
fer codes. 
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